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ABSTRACT 


A highly linear passive floating element wall shear’ trans- 
ducer has been developed which uses a sensitive frequency modula- 
tion technique to measure a capacitance change proportional to 
the displacement of the floating element for small deflections. 
The floating element is 4 mm square surrounded by a gap nominally 
0.18 mm and protrudes into the boundary layer 9 ym. Scaled on 
inner parameters, each side of the sensing head measures about 
290, the mean gap width is 13, and the element protrusion is 0./. 
The transducer was tested in a turbulent boundary layer on a flat 
plate at momentum thickness Reynolds numbers of 3140, 3500, and 
3760 and the results compared to wall shear stress estimates ob- 
tained from the law of the wall. The transducer was found to 
predict the steady mean wall shear stress within the limits of 
experimental uncertainty in all cases; transducer measurements 
of mean wall shear stress for momentum thickness Reynolds numbers 
of 3140 and 3500 were within 2 percent of the actual value. 
Transducer measurement of the unsteady mean wall shear stress was 


severely limited by a mechanical resonance in the transducer at 
OG HZ. 


The operating principles and general characteristics of the 
transducer are discussed, its performance is evaluated, and re- 
commendations are made concerning future work. 
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NOMENCLATURE 


surface area of the floating element; 
capacitance of the cantilever-stationary plate pair 
friction coefficient 


mean spacing between stationary and cantilever plates 
with the cantilever undeflected 


uncertainty or error in the quantity ¢ 
the flexural rigidity of the cantilever beam 


generally, a force; in Appendix C, a functional ex- 
pression for the quantity being measured 


frequency in Hz; in Appendix A, a force distributed 
the length of the beam element (f(x,t)) 


center frequency of astable multivibrator 


fundamental mechanical resonant frequency of the beam 
element 


gain in the frequency~to-~voltage converter and instru- 
mentation amplifier, V/Hz 


mass moment of inertia of the floating element 
surface roughness 


static spring constant of cantilever beam to transverse 
loads, 3EI/L° 





V(t) 


length of the cantilever plate and beam element 


mass of the floating element (assumed to be a point 
mass) 


mass of the’ beam 


the total head, Ne ; also the total error (Appen- 
dix C) 2 i 


electrical resistance 


Reynolds number based on x, the distance from the 
plate’ s leading edge 


momentum thickness Reynolds number 


sensitivity of the transducer in either V/Pa or V/un; 
also the spacing between capacitor plates (Appendix A) 


thickness of the cantilever plate ; also time 
longitudinal fluctuating fluid velocity 


free stream velocity 





friction velocity, 


transducer offset voltage 


transducer output voltage 


width of the cantilever plate; also (Appendix A) the 
weight per unit length of the beam 


weight of the floating element 





5x 


A\*® 


a 


10 


referring to the flow, the streamwise distance from the 
plate’s leading edge; referring to the beam element, 
the distance along the beam from the base of the beam 


referring to the flow, the distance normal to the wall 
measured from the wall; referring to the beam element, 
the deflection of the beam element from its undeflected 
position 


the deflection of the beam with respect to a stationary 


observer 


boundary layer thickness; also (Appendix A) the varia- 
tional operator used in the variational calculus 


displacement thickness 


transverse deflection of the free end of the 
cantilever; also the percent most probable error (Ap- 
pendix C) 


Vv 
the viscous length scale, a 
z 
beam deflection normalized on mean capacitor spacing d. 
normalized coordinate along beam, x/L 
momentum thickness 
transducer inclination angle from the vertical 


floating element angle of rotation 


wall shear stress; mean wall shear stress refers to a 
spatial average as opposed to steady wall shear stress 
which implies a time average 


kinematic viscosity 





+ 


Jd) 


Superscripts 


identifies the quantity as being nondimensionalized 
inner coordinates; velocities are normalized on U 
lengths are normalized on i* ‘ 


Subscripts 
refers to free stream properties 


refers to quantities measured at the wall 


in 
and 
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Ll. INTRODUCTION 


1.1 BACKGROUND 


For many centuries, men have studied the force of drag which 
resists the motion of a body through a fluid. Todd (1967) re- 
ports that as early as the late fifteenth century, Leonardo da 
Vinci conducted tests on three ship models with different hull 


forms, presumably to identify the most efficient design. 


The dawn of the “age of steam” in the late nineteenth 
century motivated researchers to better understand fluid drag in 
order to both predict and minimize its effect. As sailing ships 
were supplanted by larger and faster self-propelled ships, fric— 
tion drag was recognized as a sizeable contributor to the overall 
fluid drag. William Froude (1872) conducted a series of friction 
drag measurements in his towing tank at Torquay, England, to ana- 
lyze the effect of friction drag on ships. Froude built planes 
from wooden planks and towed them edge-first through the water, 
measuring their drag with a dynamometer. In so doing, Froude acy 
tually determined the mean wall shear stress acting over the sur- 


face of each plane. 


Since 1904, when Prandtl presented his boundary layer theory 


to the Mathematical Congress in Heidelberg, understanding of the 





wall shear stress has grown considerably in sophistication. 
While Froude intuitively grasped that the distribution of wall 
shear stress over a body in a fluid flow was not uniform, theory 
now predicts that the wall shear stress is a function of both 
space and time. This dynamic, nonhomogeneous nature of the wall 
shear stress is particularly evident in the transition and turbu- 


lent regions of the boundary layer. 


The wall shear stress is also an important boundary layer 
parameter in its own right since it is a characteristic of the 
inner region of the boundary layer (near the wall) where viscous 
forces predominate. Many researchers have found that boundary 
layer processes scale on the wall shear stress, particularly 
those related to turbulence in the transition and fully turbulent 
regions of the boundary layer. One example is the familiar ‘law 
of the wall” in which the turbulent mean velocity profile is non- 
dimensionalized on mean wall shear stress and kinematic viscosity 
[see Coles (1968) and White (1974)]. Burton (1974) found that 
turbulent wall pressure spectra and corresponding fluctuating 
longitudinal velocity spectra scaled on the wall shear stress. 
Gedney (1979) was unable to collapse his transition wall pressure 
spectra because of inaccurate determination of the wall shear 
stress upon which they were nondimensionalized. In short, ongo- 
ing fundamental boundary layer research indicates that more pre- 
cise measures of the wall shear stress than currently available 


are increasingly necessary. 





1.2 WALL SHEAR STRESS 


Wall shear stress, like any stress, is a derived quantity 
and thus cannot be measured directly. The wall shear stress is a 


dynamic quantity defined at a point by the following limit: 


oe i i on 


"yw 7 AF0 A ce) 


where F is the friction force acting over area A. Unless’ the 
measurement of ° the ea shear force is made at a single point, 
the best we can expect to measure is a mean wall shear’ stress 
averaged over some measurement area. Measurement techniques used 
to determine wall shear stress are, in reality, merely techniques 
which measure quantities related to the wall shear stress from 
which the wall shear stress is inferred. As the measurement area 
approaches a point, the measured wall shear stress approaches the 
true wall shear stress according to equation 1.1. Small measure- 
ment areas thus strengthen the inference that the measured shear 
stress is the true shear stress, underscoring the need for good 
spatial resolution if the true shear stress at a point is the 


subject of investigation. 
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1.3 MEASUREMENT TECHNIQUES 


Table I, adapted from Winter (19//), provides a brief list- 
ing of the mean wall shear stress measurement techniques current- 
ly in use. No eee technique is universally acceptable because 
each has restrictions on its application to different flow condi- 
tions and possesses limited versatility. Further, no _ technique 
can make an adequate measure of the true wall shear stress, owing 
to excessive spatial averaging or inadequate bandwidth. The hot 
film probe, using the heat transfer analogy, is the only device 
with sufficient bandwidth to measure some of the unsteady char- 
acter of the mean wall shear stress. However, its spatial reso- 


lution and sensitivity are both poor. 


TABLE I 
Mean Wall Shear Stress Measurement Techniques 


Analogies . Flow About Obstacles 
Heat Transfer Preston Tubes 
Mass Transfer Stanton Tubes 


Razor Blades 
Steps & Fences 


Liquid Films Static Pressure Holes 


Floating Element Momentum Balance 


Velocity Profiles (Wall Similarity) 





In an effort to improve the resolution of dynamic wall shear 
stress measurements, a new type of floating element transducer 
has been developed and tested. This device uses a 
capacitance-measuring technique of exceptional sensitivity [first 
developed by Moeller (1983)] to measure the deflection of a small 
floating element in the plane of the wall. The sensitivity of 
the transducer permits its operation with such small motions of 
the floating element that the device exhibits excellent lineari- 
ty. Thus, the powerful techniques of linear systems analysis may 
be used to uniquely and simply relate statistics of the transduc" 
er output to corresponding statistics of the mean wall shear 


stress. 





2. TRANSDUCER DESIGN 


2.1 GENERAL DESCRIPTION 


The transducer is a passive floating element instrument 
which generates a voltage proportional to the displacement of the 
element and, indirectly, to the shear force acting on the ele- 
ment. This is unlike most floating element devices, which gener- 
ally employ an active force-balance principle so that the element 
remains nearly stationary. The work of Kistler [Paros, (19/0)], 
Allen (1980), Frei and Thomann (1980), Escudier, et. al. (1982), 
and Nguyen, Dickenson, et. al. (1984) are all examples of the 
active technique. In the current design, sensitivity has been 
improved to the point that the floating element deflects less 
than 5 Um. For a streamwise deflection of 5 ym, the floating 
element protrudes only 0.63 ym into the boundary layer, which is 
a y+ of Belo - The protrusion due to floating element deflec- 
tion is thus well within the viscous sublayer and can be neglect- 


ed. 


Figure 2.1 presents a schematic of the transducer. The beam 
element is the moving element of the transducer, composed of the 
floating element and a thin cantilever plate. The floating ele- 
ment, mounted in the plane of the wall atop the cantilever plate, 


is sensitive to motion in only one direction. Shear forces on 





the surface of the element in the sensitive direction deflect the 
cantilever plate, producing a change in capacitance with respect 
to a second plate parallel to and mounted a short distance from 
the cantilever plate. This change in capacitance, proportional 
to floating element deflection for small displacements, changes 
the frequency of a square-wave oscillator. The transducer’ thus 
produces a frequency-modulated signal proportional to the float- 
ing element displacement which is subsequently demodulated and 


filtered by signal conditioning electronics. 


FLOW 
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FIGURE 2.1 Transducer Schematic — General 
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2.1.1 Transducer Mechanical 


The transducer is assembled within a body of 
phenolic-impregnated linen as shown in Figure 2.2. A shallow 
slot, milled in the body halves before assembly, contains a_ sta- 


tionary capacitor plate of copper-clad laminate epoxied in place. 


FLOATING ELEMENT ACCELEROMETER 


(4mm x 4mm) CANTILEVER | _ LOCATION 
(145mm x 30mm) 


Uy 


ow 





See May 
AXIS 





STATIONARY LINEN PHENOLIC 
OO PLATE BODY (2 HALVES) 
Lem (145 mm x 30 mm) 


FIGURE 2.2 Transducer Schematic 


The movable capacitor plate, made from 0.005 inch (0.13 mm) thick 
sheet brass’ shim stock, is sandwiched between the halves of the 
transducer body which are epoxied and screwed together. A float- 


ing element 4 mm square is cast in place atop the cantilever 
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plate using an epoxy potting compound. The square shape of the 
floating element was selected to minimize the degradation in 
transducer sensitivity caused by sealing the gap around the 


element. 


The transducer is then mated to an aluminum adapter machined 
to accomodate the test fixture as indicated in Figure 2.3. An 
aluminum cylinder also mates with the adapter, shielding the 
transducer and its electronics from electromagnetic interference. 
Brass shims between the transducer and the adapter are used to 
position the floating element parallel to the top of the adapter 
with the top of the floating element protruding less than 0.1 mm 
from the surface of the adapter. The tops of the adapter and 
floating element, spray painted with flat white enamel before as- 
sembly, were lightly sanded with fine wet sandpaper to reduce the 
misalignment of the floating element with the wall to 9 um. By 
sanding perpendicular to the transducer’s sensitive direction, 
motion of the floating element and curvature of its upper surface 
were minimized. The gap around the element is not symmetric, 
ranging from 0.1 mm to 0.25 mm with a mean of 0.18 mm. An accel- 
erometer mount is provided on the side of the transducer assembly 
inside the cylinder in order to monitor the acceleration of the 


transducer in its sensitive direction. 


Table II summarizes the key lengths of the floating element 


in both dimensional form and in inner coordinates, nondimension- 








Zi 


alized on the viscous length A* characteristic of the boundary 


layer in which the transducer was ultimately tested. 
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FIGURE 2.3 Transducer - Pictorial 





TABLE II 
Floating Element Dimensions 







a a 
To.10 wm [7 
To.18 mm | 
rroerusion | ova | 0.69 


Size 











gap min. 






ae 
¢ is based on the experimental result 
UL =i 3 omy Ss) at U_ = 30.5 m/s. 











2.1.2 Transducer Electronics 


The electronics associated with the transducer fall into two 
categories: electronics physically mounted beneath the phenolic 
transducer assembly, and signal conditioning electronics in- 
stalled in a separate enclosure. Each set of electronics uses 
shielded precision power supplies (Analog Devices model 904) 
which are bypassed to ground through 1.0 UF electrolytic capaci- 
tors. Further, each integrated circuit power supply pin is by- 


passed to ground through a 0.1 uF ceramic capacitor. 


Figure 2.4 highlights the significant features of the elec- 


tronic circuitry. 
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FIGURE 2.4 Circuitry - Block Diagram 





Transducer-mounted electronics consist of the frequency 
modulator anda buffer amplifier. The frequency modulator is an 
astable multivibrator which produces a square wave nominally 
6.2 V in amplitude at a frequency determined by the variable ca- 
pacitor formed by the fixed and cantilever plates. When the can- 
tilever plate is undeflected, the capacitance measures about 
30 pF and the oscillator’s frequency is about 55.6 kHz. An LM310 
voltage follower is used as an isolation amplifier that drives a 
shielded BNC transmission line to the signal conditioning elec- 
tronics. Power is supplied to the transducer-mounted electronics 
through a shielded cable from a power supply in a separate enclo- 


sure. The multivibrator circuit is taken from Jung (1980). 


The signal conditioning electronics are mounted in a shield- 
ed aluminum box containing an integral power supply. An Analog 
Devices 453L precision frequency-to-voltage converter demodulates 
the transducer’s frequency modulated signal, converting it to a 
voltage that is proportional to the received frequency. Both the 
gain and the offset of the frequency-to-voltage conversion are 
adjustable for calibration. An AD524 instrumentation amplifier 
(with selectable gains of 1, 10, 100, or 1000) amplifies the de- 
modulated signal and drives both an active filter and a transmis-— 
sion line for measuring the fluctuating signal. The active 
filter, taken from Lancaster (1982), is a first-order low pass 
filter with cutoff frequency of 0.2 Hz. The filter drives a sep- 


arate line for measuring the steady transducer output. 





ed 


2.2 TRANSDUCER RESPONSE CHARACTERISTICS 


It is shown in Appendix A.3 that the response of the trans- 
ducer is described by the following equation in the limit of 


small deflections of the floating element: 


o> Aw G £ 
Oe t_(€) Ca) 
8 Ke d. 
where 
ene) = the mean wall shear stress acting on the floating 
element; 
Vit) = the output voltage from the signal conditioning 
electronics; 
A, = the surface area of the floating element exposed 
to the flow; 
d. = the mean separation between the undeflected capa- 
citor plates; 
E = the oscillator center frequency; 
G = the total gain in the frequency-to-voltage con- 
verter and the instrumentation amplifier; 
K. = the transverse stiffness of the cantilever beam 
when statically loaded at its free end; and 
Ve = the offset voltage. 


Equation 2.1 is applicable for small displacements of the float- 





ing element and for frequencies much smaller than the first reso- 


nant frequency of the beam element. 


Equation 2.1 may be rewritten as 


VG&aas= i “+ S_ T(t) : (292) 
where 
pie, 2 “0 0.40 v/ 
> ets rs coco’ : Pa a 
T t 3 Kk A (2:43) 


is the transducer’s sensitivity. A related sensitivity is the 


transducer output per micron of floating element deflection 





Ee A as a = 0.56 V/um ° (2.4) 


The limits of the validity of equation 2.1 may be investi- 
gated by using the predicted beam element deflection shape deter- 
mined in Appendix A.3 to evaluate the theoretical transducer out- 


put. From the appendix, 


V(t) - ve : ( 235) 
ee l Seal Garis) 
G . dé 
fi renee) 
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and 





1 Gcl8) ye (la (2.6) 
O 
(A3.4) 


V(t) - V 
—°—. -- -2 the normalized transducer output andn(é,t)is 


the deflection shape of the beam element normalized on do: 
Equation 2.6 describes the static deflection shape of the beam 
element for a force applied to the floating element and is valid 
for frequencies much smaller than the fundamental mechanical re- 
sonance. The integral in equation 2.5 may be performed numeri- 
cally using equation 2.6 and the normalized transducer output 


plotted as a function of the normalized end deflection of the 


€ 
beam Fl oo 8 shown in Figure 2.5, the transducer is linear to 
O 
€ 
within 1.25 percent proportional error up to |g ;} = 0.05. Given 
O 








2 
that k, = 22.5 N/n, A, = 16 mm , and de *z 375 um, the device 


is effectively linear up to 


=| = 26 Pa since 
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For the turbulent boundary layers considered in the current’ re- 


ae ~ 35 Pa so that the normal operating regime is the 


+ 


search, - 








much smaller box shown in Figure 2.5. 


The natural frequencies of the transducer response are 
obtained in Appendix B by evaluating the eigenvalues of the beam 
element’s free vibration response. For the first two vibration 


modes of the beam element, the rotational inertia of the floating 
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element was calculated to affect the natural frequencies by less 
than 1.0 percent and so was neglected. Calculation of the trans- 
ducer’s natural frequencies is complicated by an ambiguity in the 
beam length introduced by the spatial extent of the floating ele- 
ment. Contrary to the model, the floating element is not a point 
mass but rather extends over 6 percent of the cantilever plate’s 
length. The fundamental resonant frequency of the beam element 
is estimated from the model to be 61./ Hz, with an uncertainty of 
6.3 percent. Measurements of the dynamic characteristics of the 
transducer presented in greater detail in section 4.1 indicate 


that the fundamental resonance is actually 63 Hz. 


Characteristics of the transducer are summarized in Table 


IIl. 


TABLE IIL 
Transducer Characteristics 


Mass: Floating Element 76 mg 
Cantilever Plate 480 meg 
= 16 g/m 
Stiffness: 22.5 N/m 
Fundamental Resonance: 63 Hz 
Sensitivity: 3. 0.40 V/Pa 
s = 0.56 V/um end 


deflection 


Spatial Resolution: 4 mm x 4 mm 





Inclination of the transducer from the vertical in a gravity 
field generates a response from the transducer, since then a com- 
ponent of the weight of the beam element is transverse to the 
bean. The transducer exhibits buckling behavior when inclined, 
because the transverse force tends to increase the deflection of 
the beam and itself increases as the local inclination angle of 
the beam increases. Transducer inclination response is analyzed 
in Appendix A.4 in the limit of small static inclinations of the 
transducer base. The expression obtained in the appendix for the 
deflection of the beam element n(&,t) with inclination angle 
8 (equation A4./7) is then inserted in equation 2.5, which is in 
turn integrated numerically to obtain the transducer output as a 


function of inclination angle plotted in Figure 2.6. 


The theoretical response of the transducer to inclination is 
sensitive to the beam length L and mean capacitor spacing d_ , as 
indicated in Figure 2.6 by comparing the solid and dashed lines. 
The solid line is the theoretical transducer response for the 
beam length and mean capacitor spacing used in all calculations; 
the dashed line is’ the theoretical response for a beam length 
1 mm (3.24) smaller and mean capacitor spacing 25 um (6.72) 
larger. The variations in the parameters used to generate the 
dashed line are representative of uncertainties in those parame- 


ters due to manufacturing limitations. 


A simple experiment was conducted to verify the theoretical 





Jt 


response predicted in Figure 2.6. Experimental results are also 
plotted in Figure 2.6 and may be seen to lie within the limits of 
uncertainty in predicted response. Also indicated on the figure 


are the expected operating limits of the transducer in the cur- 


rent research from Figure 2.5, corresponding to oe ae 5.0 Pa at 








zero inclination angle. The transducer output due solely to ro- 
tation of the transducer base one degree from the vertical is of 
the same order as the transducer output due solely to wall shear 
stress, indicating the importance of maintaining the attitude of 


the transducer. 
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3. EXPERIMENTAL EQUIPMENT AND INSTRUMENTATION 


3.1 EXPERIMENTAL FACILITY 


The transducer was tested in the MIT low-noise, 
low-turbulence open circuit wind tunnel. This facility is shown 
in Figure 3.1 and is described in detail by Hanson (1969). Lt 
has a flow speed range of 20 to 50 m/s, controlled by varying the 
speed of a DC motor driving the blower. The array of turbulence 
reducing screens between the honeycomb flow straightener and the 
settling zone maintains the free stream turbulence level less 
than approximately 0.054. A test section 38 cm (15 inches) 
square and 183 cm (6 feet) long was mounted in an open jet confi- 
guration within the airtight blockhouse immediately behind the 
contraction. The test section, lined with smooth sound absorp- 
tive material, has movable walls for controlling the longitudinal 
pressure gradient. A mechanical x-y traverse accurate to 0.025 
mm (0.001 inch) vertically and to 2.5 mm (0.1 inch) horizontally 
is installed on the test section and used to position a hot wire 
or a pressure probe during the measurement of mean flow proper- 


ties in the test section. 
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3.2 TEST FIXTURE 


The test fixture is a damped aluminium flat plate described 
in detail by Gedney (1979), who used it while investigating tran- 
sition region wall pressure fluctuations. The plate, sketched in 
Figure 3.2, is 38 cm (15 inches) wide and 122 cm (48 inches) 
long; the upper surface is polished to a 0.1 micron root’ mean 
square (rms) surface roughness. The plate’s leading edge is a 
solid tool and jig plate, 1.27 cm (0.50 inch) thick, machined 
into a 6:1 ellipse. The test plate’s trailing edge consists of a 
thin splitter plate which extends rearward to prevent coherent 


vortex shedding from creating plate vibrations. 


The test plate was mounted horizontally in the test section 
13 cm (5 inches) above the bottom to minimize secondary flow ef- 
fects observed by Shapiro (1977). The walls of the test section 
were then adjusted to provide a minimum pressure gradient. The 
variation in pressure obtained over the length of the plate at a 
free stream velocity of 35 m/s is shown in Figure 3.3. The maxi- 
mum variation observed (excluding the leading edge) was less than 
2% of the total head. The transducer was installed at the loca- 


tion indicated. 
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FIGURE 3.3 Longitudinal Pressure Variation 
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A flow visualization to check the effect of transverse wall 
turbulence was performed using Shapiro's technique, also used by 
Gedney. The results were identical to Gedney s as shown in Fig- 
ure 3.4. The shear stress transducer was tested in the 7.6 cm (3 


inch) hole at the location indicated. 


The boundary layer was tripped to initiate turbulence by two 
thicknesses of duct tape 1 cm in width and 0.5 m in total 
height. The tape was placed on the upper surface of the plate 
2cm from and parallel to the leading edge. According to the 
criterion developed by Gibbings [White (1976)], fully effective 


tripping to turbulent flow is assured if 


8 





= 826 (3.1) 


where k is the height of the surface roughness. In the present 
work, k =~ O.5 mm ensures immediate transition for free stream 


velocity greater than 25 m/s. 
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3.3 FLOATING ELEMENT MISALIGNMENT 


Allen (1977) has indicated that the major sources of error 
in floating element shear stress measuring devices are lip forces 
and normal forces which contaminate the friction force that the 
device is being used to measure. Lip forces are developed by a 
nonuniform pressure distribution around the periphery of the 
floating element caused by either flow through the gap around the 
element or misalignment of the floating element with respect to 
the wall. Normal forces introduce error when a nonuniform pres- 
sure distribution over the surface of the floating element shifts 
the normal force's center of effort, generating a moment which 


deflects the floating element. 


Allen found that errors in the measured wall shear’ stress 
due to lip and normal forces can be larger than the actual wall 
shear stress by as much as a factor of 4. His results indicate 
that proper alignment of the floating element in the wall is the 
key factor in minimizing error, a finding corroborated by Burke 
(1984) who used a floating element to measure skin friction in 


the MIT water tunnel. 


Proper alignment of the floating element in the wall is 
therefore critical in the current design, and was assured by 


measurement with a KD-45 Fotonic Sensor made by Mechanical Tech- 
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nology Incorporated. The Fotonic Sensor is a non-contact measur- 
ing device that projects several beams of light through fiberop- 
tic strands and then measures the quantity of light reflected 
from an object to determine the distance to the object. The Fo- 
tonic Sensor was used in its less sensitive range ( +2 um) to de- 
termine that the floating element protruded 9 ym from the wall 


(y" = 0.69), which was judged to be satisfactory. 





3.4 INCLINATION RESPONSE 


Experimental verification of the transducer’s response to 
inclinations from the vertical were made, as shown in Figure 2.6. 
The transducer was attached to the mid-~span of a sturdy 1.83 m 
long aluminum I-beam which was raised or lowered at one end a 
known distance to generate inclination angles smaller than one 
degree. An inclination angle of one degree is obtained when one 
end of the I-beam is 32 mm higher than the other end, so it is 
possible to maintain accuracy to about 0.01 degree. The effect 
of I-beam flexure on inclination angle is small to begin with and 
further reduced by placing the transducer at the I-beam’s 


mid-span, so that I-beam flexural effects have been neglected. 
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3.5 INSTRUMENTATION 


3.5.1 Dynamic Characteristics 


Transducer dynamic characteristics were investigated using 
the instrumentation shown in Figure 3.5. The transducer was 
mounted in a phenolic bracket and secured to the vibration ex- 
citer so that the transducer could be shaken either in the plane 
of the wall or normal to the wall. As shown in Appendix A.2, the 
transducer’s response to acceleration of its base and its re- 
sponse to a shear force acting on the floating element are dynam- 
ically equivalent, although caution must be exercised in inter- 
preting results since the acceleration of the transducer’s base 
also introduces an effective force along the length of the canti- 
lever plate which is absent during normal operation of the trans- 
ducer. An accelerometer, calibrated using a General Radio model 
155/7-A vibration calibrator, was attached to either the phenolic 
bracket or to the transducer body to measure the vibration of the 
transducer body. A Hewlett-Packard 5423A structural dynamics an- 
alyzer was used to generate the random noise excitation of the 
vibration exciter and also to calculate the transfer function 
between the transducer output and transducer body acceleration. 
The spectrum of the random noise excitation was shaped using an 
MXR 3l-band equalizer so that the spectrum of the transducer body 


acceleration was approximately uniform over the frequency band of 
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FIGURE 3.5 Instrumentation - Dynamic Characteristics 
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interest. 


3.5.2 Test Conditions 


The free stream velocity into the test section was monitored 
by measuring the pressure drop across the wind tunnel contraction 
section with a Betz micromanometer. A calibration was performed 


by Shapiro (1977) to ensure an accuracy of 0.1 m/s. 


The longitudinal pressure variation along the test plate was 
measured using a small static pressure tube (United Sensor and 
Control Corp. PSA-16) and a Validyne DP15-24 differential pres- 
sure transducer. A Validyne model CD23-1071 electronic manometer 
powered the pressure transducer and provided a digital readout of 


the differential pressure in millimetres of water. 


Mean velocity profile measurements were performed using a 
hot wire anemometer probe (Thermo-Systems Inc. probe model 
1261-T1.5) with a Disa model 55005 battery-powered constant tem- 
perature anemometer and a linearizer built by Moeller using an 
Analog Devices 433B programmable function generator. The hot 
wire probe was calibrated against a pitot tube used with the Val- 
idyne DP15-24 differential pressure cell. Mean velocity was 
measured ona digital voltmeter after being filtered by a 0.1 Hz 


active low pass filter built by Lueptow (1984). 





46 


3.5.3 Transducer Calibration 


Wall shear stress transducer electronics were calibrated 
using a General Radio model 1210-C function generator in conjunc— 
tion with a Hewlett-Packard model 5532A electronic counter and a 
digital voltmeter. With gain of the AD524 instrumentation an- 
plifier set at 10, the gain of the frequency to voltage converter 
was adjusted so that a change in input frequency of 500 Hz (from 
55.6 kHz to either 55.1 kHz or 56.1 kHz) produced a change in 
output voltage of approximately 5.0 V. The exceptional sensitiv- 
ity of the transducer made workbench calibration of the entire 
transducer unfeasible, so the transducer was calibrated against 


the law of the wall in situ. 


3.5.4 Wall Shear Stress 


The wall shear stress was measured using the instrumentation 
of Figure 3.6. Coherence between the acceleration of the trans- 
ducer and transducer output was monitored using the structural 
dynamics analyzer to ensure that transducer acceleration did not 
significantly affect transducer output. The structural dynamics 
analyzer was also used in its time record mode to make 


long-period zero drift measurements. 
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4. EXPERIMENTAL RESULTS AND DISCUSSION 


4.1 DYNAMIC CHARACTERISTICS 


Transducer dynamic characteristics were evaluated for two 
different transducer configurations: one with a seal in the gap 
surrounding the floating element and another with the gap 
unsealed. For brevity, these two configurations will be referred 
to as the “sealed” and “unsealed” configurations, respectively. 
Investigation of the dynamic characteristics of the two gap con- 
figurations was conducted in two phases. In the first phase, the 
accelerometer used to measure transducer body acceleration was 
attached to the phenolic bracket mounting the transducer to the 
vibration exciter. This test configuration permitted the rota- 
tion of the transducer in the bracket and was used to verify that 
the transducer can discriminate against signals orthogonal to its 
sensitive axis in the plane of the wall, but introduced errors in 
that the accelerometer was not mounted directly on the transducer 
body. The second phase utilized an accelerometer attached di- 
rectly to the transducer body in the location indicated on Figure 
2.2. In this later phase, attention was devoted exclusively to 
improving understanding of the dynamics of the transducer in its 
sensitive direction, the orthogonality of the device having 


already been demonstrated. 
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Phase 1. The sealed configuration in the initial phase used 
as the gap sealing material an epoxy, ECCOGEL 1365-15, made by 
Emerson and Cuming. Measuring only 15 on the Shore’ hardness 
scale commonly used in the polymer industry, the epoxy is very 
soft but is nearly incompressible. Owing to the manner of the 
transducer’ s construction, it was not possible to determine the 
geometry of the gap seal beneath the surfaces of the adapter and 


floating element after the epoxy cured. 


The transfer function between the transducer output and the 
transverse acceleration of the transducer body is shown in Figure 
4.1 (unsealed configuration) and Figure 4.2 (sealed configura- 
tion). The magnitude plots are normalized on the peak magnitude 
of the unsealed transducer’s transfer function in order to high- 
light the order of magnitude loss in response caused by the epoxy 
seal in the gap. The fundamental resonant frequency of the 
unsealed transducer (Figure 4.1) is approximately 63 Hz, very 
close to the 61./7 Hz predicted during design. Figure 4.2 indi- 
cates that the addition of the epoxy seal leads to an increase in 
the fundamental resonant frequency to about 360 Hz as well as an 


apparent reduction in sensitivity. 


The apparent loss of sensitivity is misleading. Table V 
(Appendix B) lists the resonant frequencies for the first three 
vibration modes of the beam element subjected to clamped-free, 


clamped-pinned, and clamped-clamped boundary conditions. The ob- 
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served fundamental resonant frequency of the sealed transducer 
(360 Hz) is bracketed by that of a clamped-pinned beam element 
(341 Hz) and a clamped-clamped beam element (494 Hz), neglecting 
the rotational inertia of the floating element which was found to 
be insignificant. This suggests that the beam element in the 
sealed configuration exhibits a pinned condition at the floating 
element or, more correctly, a mixed condition that is a hybrid of 
pinned and clamped conditions but closer to pinned. If such is 
the case, then the sealed transducer does not respond to wall 
shear stress, since the floating element is located at a node of 
the beam element”’s response and therefore never translates. The 
sealed transducer is thus incapable of responding to forces ap- 
plied to the floating element, but rather responds only to forces 
applied to the cantilever plate. This is shown in Figure 4.2, 
where the transducer is responding to the D°“Alembert or inertia 


d¢z 


force  _ (using the notation of Appendix A.2) produced on 


us” 
dt 

the cantilever plate by the transverse acceleration of the trans- 

ducer base. The use of a seal in the gap around the floating 


element is self-defeating in this instance, eliminating transduc- 


er response to the wall shear stress. 


The resonances of both the sealed and unsealed configura- 
tions exhibit considerable damping, approximately 14 percent of 
critical damping in each case. The unsealed configuration’s 
response in Figure 4.la shows an irregularity at 60 Hz caused by 


noise from the field windings of the vibration exciter, while the 





poor response below 10 Hz is due to a low frequency limitation of 
the vibration exciter. The phase angle of both configurations” 
transfer function is anomalous, since it does not change by 180 
degrees as it should when passing through an isolated resonance. 
The irregularity in Figure 4.la near 50 Hz and the similar irre- 
gularity in Figure 4.2a near 325 Hz are also strange, but may be 
due to the vibration of the bracket securing the transducer to 


the exciter. 


Phase 2. The second phase of the investigation of transduc- 
er dynamics was undertaken to resolve some of the issues raised 
concerning the phase of the transfer function and the role of the 
mounting bracket in introducing the anomalies encountered in the 
earlier phase. Also of interest was confirmation of the hypo- 
thesis that the floating element in the sealed configuration was 
effectively held motionless by the seal. In this phase, the ac- 
celerometer used to measure transducer transverse acceleration 
was mounted directly on the phenolic transducer assembly rather 
than on the bracket mounting the transducer to the vibration ex- 
citer. For the “sealed” configuration, a single piece of stan- 
dard vinyl electrical tape was securely attached to the surfaces 
of the floating element and adapter to approximate a pinned boun- 


dary condition at the floating element. 


Figures 4.3 and 4.4 illustrate the transfer function 





obtained between the transducer output and the transverse accel- 
eration of the transducer body. The units of the ordinate in the 
Magnitude plot are volts per g, where l g = 9.8 m/s” is the ac- 


celeration due to gravity. 


The unsealed configuration’s transfer function (Figure 4.3) 
is virtually identical to Figure 4.1, with the following impor- 
tant exceptions: 

- the phase is well-behaved and very properly shifts 
through 180 degrees when passing through the resonance; 
- the anomaly in the magnitude near 50 Hz has vanished, 
as expected; and 
- the irregularity in magnitude at 60 Hz is sharper’ than 
that encountered earlier, and gives rise to slight ir- 
regularities at 120 Hz and 180 Hz. 
As discussed earlier, the noise at 60 Hz (and multiples thereof) 
is due to the field windings in the vibration exciter. The 
slight irregularity in magnitude at 100 Hz is due to noise from 


the vibration exciter’s blower motor. 


The "sealed" configuration’s transfer function is remarkably 
similar to that of Figure 4.2, despite the difference between the 
two seals (electrical tape as opposed to epoxy). In Figure 4.4: 

- the fundamental resonant frequency has increased _ to 
375 Hz (from 360 Hz in Figure 4.2); 


- the damping has decreased to 7 percent of critical (in 





MAGNITUDE (V/x) 


PHASE ANGLE (DEGREES) 


120 


100 


80 


63 


0 20 40 60 80 100 120 140 160 


FREQUENCY (HZ) 


a. Magnitude 


180 


I> 


90 


~45 


00 


-135 


63 


-180 


0 20 40 60 80 100 120 140 160 


FREQUENCY (HZ) 


b. Phase Angle 


180 


180 


200 


FIGURE 4.3 Final Transfer Function - Gap Unsealed 





MAGNITUDE (V/g) 


PHASE ANGLE (DEGREES) 


> 
© 


3/5 


0 100 200 300 400 500 600 799 


FREQUENCY (HZ) 


a. Magnitude 


135 


90 


-135 
a5 





=130) 
0 100 U0 300 400 500 600 700 


FREQUENCY (HZ) 
b. Phase Angle 


FIGURE 4.4 Final Transfer Function - Gap Sealed 


800 


800 





I/ 


contrast to 14 percent in Figure 4.2); 

- the large irregularity in magnitude just below reso- 
nance has vanished; and 

- the phase is well-behaved and shifts through 180 de- 


grees through the resonance. 


The results of this phase of the investigation indicate that 
the anomalies encountered were indeed experimental and confirm 
the hypothesis that the gap seal immobilizes the floating element 


in translation. 





4.2 WALL SHEAR STRESS 


4.2.1 Law of the Wall 


The law of the wall was used to obtain estimates of the mean 
wall shear stress for comparison with the steady mean wall shear 
stress measured by the transducer. Mean velocity profiles were 
measured above the floating element with a hot wire probe at free 
stream velocities of 26.0 m/s, 30.5 m/s, and 34.5 m/s correspond- 
ing to momentum thickness Reynolds numbers of 3140, 3500, and 
3760. The velocity profiles were analyzed using the FORTRAN pro- 
gram BLPROF2 in which the user iteratively guesses the friction 
velocity ue and the location of the wall (in the measured coor- 
dinates) until the measured data coincides in the inner region of 
the profile with Musker’s explicit formula for the turbulent mean 
velocity profile [see Musker (1979)]. BLPROF2 permits the user 
to fair the measured data into the wall using Musker’s formula 
and then integrates the resulting faired profile using trapezoi- 
dal rule to obtain displacement thickness 6* and momentum thick- 


ness 6. 


The three mean velocity profiles are plotted in Figure 4.5 
nondimensionalized on inner parameters of the boundary layer. It 
can be seen that the velocity profile for Re, = 3500 contains 
little or no wake, whereas the other two profiles exhibit a no- 


ticeable wake. The velocity profile at Re, = 3500 was measured 
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FIGURE 4.5 Mean velocity profiles in semilogarithmic coordinates. 


Ordinate scale refers to lowest curve; other curves 
successively displaced upward by 5.0 units. 


4, 





during a different test period than the velocity profiles at 
a = 3140 and 3760, highlighting the difficulty of reproducing 
test conditions exactly. Since the anomaly in the wake at 
Re, = 3500 does not affect the logarithmic overlap region, the 
change in the outer flow causing the different wake is expected 
to have little effect on the wall shear stress. Because all 
three mean velocity profiles collapse to the classical turbulent 
mean velocity profile, it is assumed that each boundary layer ex- 
hibits fully developed turbulence despite being tripped into tur- 


bulence so that the application of the law of the wall is valid. 


4.2.2 Transducer Measurement 


Steady. Steady mean wall shear stress was measured by _ the 
transducer at the same Reynolds numbers as the velocity profiles; 
the results are listed in Table IV, where they are compared to 
the results obtained from the law of the wall. Both the law of 
the wall data and the transducer-measured data are plotted in 
Figure 4.6, taken from skin friction data analyzed by Murlis, 
et. al., (1982). Transducer measurements at Reg = 3140 and 3500 
are excellent, while that for Re 4 = 3/60 is a bit high but still 
within experimental uncertainty, estimated in Appendix C to be 66 
percent. The acceleration of the transducer was neglected in an- 
alyzing transducer output, since low coherence between the _ two 


(at most 10% to 20%) indicated that there was little systematic 








ol 


relationship between the transducer output and acceleration of 


the transducer body. 


TABLE IV 
Results 
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Transducer zero drift was surprisingly low despite the 
transducer’s extreme sensitivity, averaging 20% of the mean sig- 
nal over the period of the measurement (1 to 1.5 minutes). Zero 
drift was measured by covering the floating element with an index 
card while the wind tunnel was running at the proper. speed. A 
short length of wire 0.5 mm in diameter was taped beneath the 
card as a spacer to keep the card from contacting the floating 
element. Sample plots of the zero drift measured by the struc” 
tural dynamics analyzer over a period of 1 minute are shown in 
Figure 4.7 where they are compared to the mean output voltage de- 
termined by subtracting the average reference voltage from the 
transducer total output voltage. Each data run took 3 to 4 
minutes from the start of the initial zero measurement to the 
measurement of the wall shear stress and the concluding zero 


measurement. 
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FIGURE 4.6. Skin friction coefficient variation with momentum thick- 
ness Reynolds number [from Murlis, et. al., (1982) | : +, 
measured data using wall shear transducer; o, data from 
law of the wall. 
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Unsteady. The unsteady mean wall shear stress Spectrum for 


Re, = 3140 and 3760 is shown with the background spectrum in Fig- 


9 
ure 4.8. Although the plots show some dynamic character of the 
unsteady wall shear stress, they show little of its structure 
owing to the low resonant frequency of the transducer. If it is 
assumed that an instrument is usable up to a frequency one order 
of magnitude less than its first natural frequency, then the 
transducer can only be used up to 6.3 Hz. Burton (1974) indi- 


cates that the frequency related to sublayer bursting in a turbu- 


lent boundary layer is given approximately by 


Ey: 0.03 (4.1) 





where 6* is the displacement thickness. For the conditions under 
which the transducer was tested, equation 4.1 predicts the burst 
frequency to be about 300 to 500 Hz, well above the nominal 


bandwidth of the transducer and thus unmeasurable. 
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Background 
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FIGURE 4.8 Wall Shear Stress Spectral Density 
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>. RECOMMENDATIONS 


In the short term, the wall shear transducer is usable for 
measuring the steady mean wall shear stress. The sensitivity of 
the device is so great, however, that calibration must be per- 
formed in situ, preferably in a laminar boundary layer where the 
wall shear stress is well known. The electronics should be up- 
graded by the installation of precision instrumentation 
amplifiers to replace the commercial-grade amplifiers currently 
in use. This should reduce the high frequency noise observed at 
the transducer output and improve the stability of the astable 
multivibrator, whose center frequency was observed to drift from 


55.6 kHz to 46.5 kHz over seven days of continuous operation. 


Long range plans for the transducer should include increas- 
ing its bandwidth, decreasing its sensitivity to inclination, and 
gaining an understanding of how the gap around the floating ele- 
ment affects the shear stress measurement (given that the gap 


cannot be sealed, as the current research suggests). 


As regards the bandwidth, the first resonant frequency of 
the beam element needs to be increased at least two orders of 
magnitude so that the transducer can measure frequencies near the 


burst rate. A larger bandwidth is desirable so that the device 
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can always be used more than an order of magnitude below its fun- 
damental resonance, which is far less complicated than attempting 
to use an instrument through its first resonance. The transducer 
should be as_ simple to use as possible. The linearity of the 
transducer, coupled with sufficient bandwidth, will make it a 
powerful instrument for measuring the dynamics of the wall shear 
stress. As linear system theory demonstrates, only in a linear 
system can the statistics of an output (such as transducer output 
voltage) be related definitively and unambiguously with statis- 


tics of the input (such as wall shear stress). 


Next, the sensitivity of the transducer to inclinations from 
the vertical should be reduced by at least two orders of magni- 
tude so that inclination has a second-order effect on transducer 
output (it is assumed that the steady transducer signal is 
zeroth-order while the unsteady signal is first-order). Coupled 
with the need to maintain roughly the same sensitivity to wall 
shear stress $3 the requirements for an improved physical design 
imply that the mass of the floating element and capacitor plate 
Must be reduced several orders of magnitude. If the size of the 


floating element is reduced, then equation 2.3 indicates that 
fo G 

ks do 
desired S_: Decreases in Ke require further decreases in mass in 


must be increased a compensating amount to maintain the 


order to maintain the required bandwidth. Increases in G lead to 
only marginal improvements in transducer performance since noise 


is generally amplified along with the signal. 
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Finally, and most important, a better understanding of the 
role of the floating element gap in affecting the boundary layer 
must be reached. There is little in the literature that sheds 
any light on this issue. Allen (1977) identified that small gaps 
increase the sensitivity of the floating element transducer to 
errors introduced by misalignment of the floating element with 
the wall. Winter (1977) mentions that the effect of gaps perpen- 
dicular to the flow is poorly understood, but that the effect of 
the gap should be negligible if the gap width is smaller than 100 
in wall coordinates. McMichael, et. al., (1979) found that a 
transverse slot in the wall of a turbulent flow affects the boun- 
dary layer downstream of the slot. In the authors” words, “Evi- 
dently, temporary release of the no slip condition as the bounda- 
ry layer passes over the slot results in a local decrease in tur- 
bulence production as indicated by the weakened intensity near 
the wall just downstream of the slot.” (p. 417). Since the gap 
is adjacent to the wall area where the shear stress will be meas- 
ured, it is unclear whether or not the boundary layer will “heal” 
itself sufficiently to permit an unbiased measure of the wall 
shear stress, particularly if the measurement area is small. The 
effect of the gap is a fundamental issue which has not yet been 


adequately addressed. 


The most fruitful course of action appears to be the micro- 


miniaturization of the transducer to a scale where a very small 
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gap and precise alignment of the floating element can be mainta- 
ined. Peterson (1982) discusses the merits of unicrystalline 
silicon dioxide at great length; small cantilever beams 8 im 
long, 2 um wide, and 0.1 um thick with a fundamental resonant 
frequency of 1.2 MHz have already been made from the material. A 
microminiaturized wall shear stress transducer with such a small 
cantilever, using the capcitance-measuring technique outlined in 
the current research, has a number of distinct advantages. A 
floating element with sensing area 2 um by 0.1 um would have su- 
perb spatial resolution, excellent sensitivity, and a large 
bandwidth. The considerable reduction in mass would greatly re- 
duce the sensitivity of the device to inclination. Further, it 
should be possible to keep the gap around the floating element on 
the order of about 10 angstroms so that the gap has a negligible 
effect on the flow. Precise alignment of the floating element is 
possible on the microscopic scale, which may circumvent the ex- 
treme sensitivity of the floating element to misalignment~-induced 
errors when the gap is small. Recent trends suggest that such a 
transducer could be mounted on a 20 um chip with its own acceler- 
ometer and compensation network, permitting an array of such 
transducers to measure the spatial distribution of the shear 


stress over the surface of the wall. 





6. CONCLUSIONS 


The development of a passive floating element wall shear 
stress transducer using a sensitive capacitance-measuring 
technique has been presented. In particular, the viability of 
using a frequency modulation scheme to perform the capacitance 
measurement has been demonstrated. The electronic sensitivity is 
such that the transducer is able to measure a 0.1 um displacement 
of the floating element. Expected motions of the floating ele- 
ment are less than 5 im and should not adversely affect the char- 
acter of the boundary layer in which the measurement is made. 
The device is highly linear, enhancing its attractiveness for dy- 


namic wall shear stress measurements. 


The transducer has been tested and found to predict the 
Steady mean wall shear stress to within the limits of 
experimental uncertainty (6642); two of the three data points ob- 
tained were within 2 percent of the actual steady mean wall shear 
stress. The transducer must be calibrated in situ in order to 
reduce the experimental uncertainty; workbench calibration is 
not feasible because of the transducer’s exceptional sensitivity. 
The current transducer is not practical for dynamic measurements 
of the mean wall shear stress since it does not possess adequate 


bandwidth. 





fs 


The gap around the floating element has been verified as 
being the greatest single vulnerability of the floating element 
transducer. It is not possible to seal the gap, since a gap seal 
in the current design immobilizes the floating element and pre- 
vents transducer response to wall shear. stress. On the other 
hand, an open gap affects the boundary layer and introduces er- 
rors of its own to the wall shear stress measurement. 
Fundamental research is needed in this area in order to fully 


understand the gap effect. 


Finally, a course of action has been recommended to improve 
the capabilities of the wall shear transducer so that it is prac” 
tical for dynamic measures of the shear stress. In order to in- 
crease bandwidth, reduce transducer sensitivity to inclination, 
and minimize the adverse effects of the gap, the transducer 
should be microminiaturized. Microelectronics technology is be- 
coming increasingly capable, and can already maintain the neces- 


sary sub-micron tolerances. 
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APPENDIX A 


CALCULATIONS 


A.l Transducer Model 
A.2 Beam Element Equation of Motion 
A.3 Transducer Response 


A.4 Transducer Sensitivity to Inclination 
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A.l Transducer Model 


Using the model of Figure A.1l, it is possible to develop a 
general theory to describe the output of the wall shear transduc— 
er with arbitrary input. Figure A.1l highlights the important fa- 
cets of transducer operation: 

- the wall shear stress causes a deflection of the beam 
element; 
- deflection of the beam element alters its capacitance 


with respect to a stationary parallel plate; 





FIGURE A.l Transducer Model 
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= the change in capacitance of the beam 
element-stationary plate changes the oscillation fre- 
quency f of an astable multivibrator, since the oscil- 
lation frequency is inversely proportional to the capa- 
citance; and 

- the frequency difference between the oscillator signal 
and a reference frequency E eof is then converted to a 
voltage which is added to an offset voltage ve » gener- 


ating the transducer output. 


The capacitance of a parallel plate capacitor is given by 


we L 
Ss 


C=K [Hudson (1944)| (Al.1) 
where K is a constant of proportionality, L is the length and w 
the width of each plate, and s is the distance between plates. 
If one plate is deflected a small transverse amount, then _ the 
distance s is no longer constant but rather a function of x, the 
coordinate along the undeflected length of the movable plate 
measured from the wall. Further, the deflection may also be a 
function of time. Effectively, s not a constant leads to an in- 


finite number of parallel capacitances 


Ca. beg? 
C=K GRD) dx ( ) 


in the limit of small plate deflections. If now s(x,t) is de- 
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fined as 


s(x,t) = d+ Veit). 5 (Al.3) 


where dis the undeflected capacitor plate spacing and y(x,t) is 


the transverse deflection of the movable plate, then 








Al 
W = w_L L 
dC = K ae y (x, t) dx = K A (x t) (A1.4) 
O o lt A 
O 
Nondimensionalizing and defining 
eae 2 = n(é,t) (Al. 5) 
L d 
O 
and noting that 
Wa 2 
K a C (Al. 6) 


where C. is the capacitance of the stationary-movable plate pair 
with the movable plate undeflected, equation Al.4 becomes 
dé 
ac = C ———— 
ol+#né,t) ° ao 
Equation Al.7 is valid only for small deflections n(&,t), for 
which three dimensional effects are insignificant. The total ca~ 


pacitance may be found to be 
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l 
: z ws ae 
C= [« = iE Pa ERE) (A1.8) 
0 


“summing” over all the infinitesmal capacitances dC. 


The oscillation frequency of the astable multivibrator is 


given by 


1 
—= ny Be 


where C is the capacitance of the transducer and Ris a_ feedback 
resistance (constant). Then, for a deflection of the movable 


plate n(&,t), the oscillator frequency is 





i ad ae 1 il (Al.10) 


RC RC ae 
0 1 + n(é,t) 


or 


ae 
7 1 
f = a ae (Adina) 
ate est) 
0 
Since 
V = _ 
(t) ee +G(f ae ) Cle?) 


it is appropriate t elect =fet btai 
ppropriate to s ect ft o¢ 5 =? 0 ain 





Vs 


1 
V(t) = ue + G CA (=) = l (Al.13) 


or 


1 
V = VY + l 
(t) Gf 


dE OE 
i + mice. &) 


using equation Al.1ll. 


The normalized transducer output is given by 


= l 
V(t) ue 


= 1 = ae (Al.15) 
Gf dé 
z 1 + n(é&,t) 


0 


The deflectionn(é,t)is determined by the transducer’s excitation 


and dynamics considered in the following section. 





A.2 Equation of Motion 


The dynamics of the movable plate of the transducer deter- 
mines how the transducer will respond to applied loads. The mov- 
able plate is cantilevered at its base and has a floating element 
mounted on its free end. The movable plate and floating element, 
together called the beam element, may be modeled as in Figure 


A.2. 


| F(t) 


=z txt) — 2( 0.) 
fe vo 


mee 





| 


“* @ z(0,t) = z(t) 


FIGURE A.2 Beam Element Dynamic Model 





The floating element is modeled as a mass M and mass' moment 
of inertia J lumped at the free end of the cantilever plate, mo- 
deled as a uniform beam of flexural rigidity EI and mass per unit 
length wy. The motion of the beam element is described by z(x,t) 
in a stationary reference frame. In order to develop the equa- 
tions of motion for general excitation, a concentrated transverse 
force F(t) is applied to the floating element and a distributed 


transverse force f(x,t) is applied along the beam. 


Of interest in the current situation is the motion of the 
beam element relative to its base rather than its absolute motion 
z(x,t), since the stationary capacitor plate is actually fixed 
with respect to the cantilever’s base. In effect, it is desired 
to obtain the equations of beam element motion in the 
non-inertial reference frame translating with the cantilever’s 


base. 


If 


z(t) =eziCOt ) (Agee) 


is the translation of the base of the beam element, then the 


relative motion of the beam with respect to its base is 





WiCceee= ZC x,t) = z(t). 


Using the variational technique [Crandall, et. 


assuming small deflections of the beam, the equation 


the beam element is obtained as 


u d*z 


92 3 oO _ 
u ae re aut ee icct) 


Q 


subject to the boundary conditions 


3 
Sy [fe ay = 0 


92 dez | 93 
by Se 1 ore + F(t) |= 0 


oy \I_ 9 Sy _ 94 _ 
8 (3) _ dt“ ox au i 0 


(A2. 2) 


and 


of motion of 


(AZ) 


(A2.4a) 


(A2.4b) 


(A2.4c) 


(A2.4d) 


where 6 is the variational operator. For the case where the beam 


is cantilevered at its base, 


Oy Set) bi) 
x=0 i x=0 


i) 
© 


so that boundary conditions A2.4a and b are identically 


fied. 


(A2.5) 


satis- 





Clamped-Free Bean. 


In the case of normal transducer operation, the beam element 


3 
is free at x=L so that the variations in y and peat x=L do not 


vanish identically, requiring the terms in brackets in equations 


A2.4c and d to be equal to zero. Thus, 


92 a d°2, (A2.6) 
ugee + BEG - fet) - uae : 
with boundary conditions 
UO Ais) = O (A? 7a) 
2 (9. £) = QO (AZ D) 
ox 
d2y(L,t) 93y(L,t) dz (A2.7c¢c) 
—3 ’ = 2, ’ 2 7 poe ne 
M JE EL a F(t) M nee 
9¢y(L,t) 93y(L,t) 
= EL = er . (A2.7d) 


Note that acceleration of the transducer base introduces 


D°Alembert or inertial force excitations proportional to 
d*z | 
O in the field equation (A2.6) and the force boundary condi- 


~ dt? | 

tion (AZ a COG If the transducer base is  non-accelerating 
d2z 

Ga = Q ), the transducer is kept vertical so that f(x,t) = 0, 


and the floating element is properly aligned so that there are no 
error forces on the floating element (F(t) = 7 OVE then the 


field equation becomes 





ait 
ust a BIL ea = 0 (A2.8) 


subject to the revised boundary condition (A2.7c) 


32 ¥(L,t) eae. = 
M = . a = 1 (t) A, (A2.9) 


with the other three boundary conditions unchanged. Therefore, 
in the absence of the error terms specified above, the wall shear 
stress excites the homogeneous response of the beam y(x,t) 
through the force boundary condition, equation A2.9. Transverse 


acceleration of the transducer base introduces not only a 


d2z 
D°Alembert force at the floating element (x=L) equal to aM ea : 
d?z 
but also a distributed D°Alembert force a along the length 


of the beam. Thus, the response of the beam element to 
transverse acceleration of its base is dynamically analogous to 
its response to the wall shear stress except that transverse ac- 


celeration gives rise to the aforementioned distributed 


D°Alembert force. 





Of 


Clamped-Pinned Beam. 
If the beam element is pinned at the floating element, then 
Sy = O at x=L and the equations describing the motion of the beam 


element are 


a 
42 a4 d ZO 
—# E = a A2. 
Lev ecuet EL ax ag a Way ( 6) 
with boundary conditions 
Cae ae (A2.7a) 
dy(0,t) - 4 (A2.7b) 
ox 
ye ee Co) 
UE aay Glee ©) 
- EI = J seo , (A2.7d) 


Because F(t) does not excite the transducer, the wall shear 
stress cannot either. Thus, if the floating element exhibits a 
pinned end condition, the beam element will not respond to the 
wall shear stress although it will respond to the tranverse acy 


celeration of the transducer base (equation A2.6 applies). 
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If the floating element exhibits a clamped end condition, 


then equation A2.6 is still the appropriate field equation, with 
y and 2¥ identically zero at x=0 and x=L. Again, the beam ele- 
ment will not respond to the wall shear stress but will respond 


to transverse acceleration of the transducer base. 
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A.3 Transducer Response 


For frequencies much smaller than the fundamental resonant 
frequency of the beam element, the deflection shape of the beam 
element is satisfactorily described by its static deflection 
shape. Assuming that the transducer base is non-accelerating, 
that the transducer is kept vertical, and that there are no error 
forces on the floating element, equation A2.8 may be solved with 
equations A2./a,b, and c and equation A2.9 for the static deflec- 


tion shape (which is not a function of time). Assuming 
Ge) A. t,t) CAs) 


leads to 


A 
e 


meet) oa ade (352 - enc) (A3.2) 
S O 


where k = 3EI/L3is the effective stiffness of the beam in re- 
S 

sisting aan applied transverse force at the free end. 

alternatively, since 


A 


c(t) = a cence) (A3.3) 
S 


is the transverse deflection of the free end of the bean, 
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ite, t) = a (ie = —-) c(t) (A3.4) 
O 


describes the deflection shape of the beam element. Since 


V(t) - Vv | 
O i dP 
———_° = Y ome 


Gf. [ dé 
(eee ACs, c) 
0 


it is possible to evaluate the transducer output as a function of 


j (Al. 15) 


T(t) 


A.3.1 Linearized Response 








Iie) Cet) | << 1, then 
1 1 1 
1 a —— a a 1+ frcese de se {A3.5) 
dE 
ea eS free dé 0 
0 
Thus, 
V(t) - Ve ; 
8 = free dé (A3. 6) 
Gf 
oO 0 


from which, with equation A3.2, 


om 1 
V(t) - V Ae Ce) 
2S 2 2 f[(3e2 = 63)ae . (A3.7) 


G f 2. Kae a 
O Ss oO 0 


Upon integration, equation A3./7 leads to 





mk 


a T(t) (A3.8) 


E 
a (A3.9) 


the normalized end deflection of the bean. 


A.3.2 Exact Response 


Using equation A3.4, the integral in equation Al.15 can be 


€ 

computed numerically for different qa: The results of this in- 
O 

tegration are shown in Figure 2.5. 





Be 


A.4 Sensitivity to Inclination 


When the transducer is inclined from the vertical by an 


angle a , the beam element deflects by its own weight. For 


small inclinations from the vertical, the transverse component 


/ la 
s,/; 
4 


Y=6 +6 
oO 






go + 


| y (x) 


beam weight per unit 
length w 


FIGURE A.3 Transducer Inclination Geometry 





fv 


i of a load F is given by 


E. er Simaves Fi) < (A4.1) 


where ¥ is the local angle between the tangent to the beam and 


the vertical shown in Figure A.3. Now, 


p= 0 +0569 + (A4.2) 


d 
where 98 ae is the local angle between the tangent to the beam 
and the undeflected beam position and os is the angle the tran- 


ducer base is inclined from the vertical. 


The beam element consists of a floating element of mass M 
and a cantilever plate of mass per unit length u. Ina gravity 
field, these possess weights of W and w, respectively. If the 
transducer is inclined from the vertical with no shear stress ap- 
plied to the floating element and without accelerating the base 
of the transducer transversely, then, from equations A2.6 and 
A2.7 and by considering only static deflections y(x) of the beam 


element, 


mr 
ee ot) (A4.3) 


x 


subject to the boundary conditions 





itl 
© 


y (0) 


(A4.4a) 
dy(0) = 4 
dx (A4.4b) 
3 
-er $3) = F(t) (A4.4c) 
dea) (A4.4d) 
ae Of 


Both f(x,t) and F(t) are developed by the weight of the canti- 


lever plate and the floating element according to equations A4.1 


and A4.2: 


- dy 
f (cet ) w ( ae + oe ) (A4.5) 
es dy 
EGE) =i Negaarpa (A4. 6) 


Equations A4.3 through A4.6 are then solved to obtain 


WS) = 74+ Gy ‘exp( ss) + 
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0 


exp(- 5 E) Lo: sin BE + C, cos oi) - & =e , (A4.7) 
O 
wL 3 
where 6 is in radians and 8° =F; ° The coefficients C, were 
calculated using the FORTRAN program GRAV3, owing to the complex- 
ity of the algebra and in order to facilitate the ready calcula- 
tion of deflection shapes for different transducer parameters. 
The transducer parameters used for the calculations are: 
L = 31 m, 


d, = 375 1m , 





2) 


M = 76 ng, 
UH = 16 g/m, and 


EL = 2.84x10 wm 


for which 
C, = 1.1568 
C = 0.9807 
2 
C = 2.1674 
5 
CG e¥62=221375 
rn 
and 8 = 0.254314, 





APPENDIX B 


TRANSDUCER NATURAL FREQUENCIES 


B.l1 Clamped-Free Beam 
The homogeneous form of equations A2.6 and A2./7 may _ be 


solved in order to determine the natural vibration frequencies of 


the clamped-free beam element: 


2 
u SF + ELSE = 0 (B1.1) 


subject to the boundary conditions 


y(0,t) = 0 (B1l.2a) 

3y(0,£) aan (B1.2b) 
2 

SY, t) Soot ets) aan (B1.2c) 
2 

Se cae) ai ge) 5 (B1.2d) 


Separating variables and assuming time harmonic behavior leads to 


yx, te —axGoae TCL) (B1.3) 


where 
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X(x) = C, sinh Ax + Cy cosh Ax + C3 sin Ax + Cy, cos Ax 


and Che a 


with 


Imposing the boundary conditions (Bl.2) obtains 


cosh AL + cos AL = y(sinh AL = sin AL) M 


sinh AL - sin AL - y(cosh AL = cos AL) iz u i 0 
where 
w* J 
sinh AL + sin AL - VET (cosh AL - cos AL) 
vy = 
wr J . 
cosh AL + cos AL - jET (Sink AL + sin AL) 


The roots of equations Bl.6 determine the eigenvalues 


(Bl.4a) 


(B1.4b) 


(Bi) 


(Bl. 6a) 


(Bl. 6b) 


i 
1 


The first three natural frequencies have been calculated and are 


tabulated in Table V. The effect of the floating element’ s 


mass 


moment of inertia was approximated by modeling the element as a 


slab 4 mm by 4 mm and 2 mm thick, and was found to be negligible 


for the first three vibration modes. 
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B.2. Clamped-Pinned Beam 
The determination of the resonant frequencies of the 


clamped~pinned beam follows the procedure for the clamped-free 


beam except that equation Bl.2c is replaced by 


Geete = a0 (B2.1) 


The characteristic equation then becomes 


Sinh AL - sin AL = y(cosh AL - cos AL) = 0 GBee2) 


where Y is as given in equation Bl.6b. The natural frequencies 
of the first three clamped-pinned modes is tabulated in Table V, 
both including and neglecting the effect of the rotational iner- 
tia of the floating element. Blevins” (1979) tabulation of the 
natural frequencies for the clamped- pinned beam are the same as 


the natural frequencies predicted from equation B2.2 above. 





For the clamped-clamped case, the mass moment of inertia no 


longer enters the problem since 


oyCent)) = 20 (B3.1) 
ox 


is the new boundary condition that replaces equation Bl.2d. The 
natural frequencies of the first three clamped-clamped modes are 


taken from Blevins (1979) and shown in Table V. 





LUU 


TABLE V 


summary of Beam Element Natural Frequencies 


Clamped-Pinned Beam f ., 


i! 5 


340.6 1104 2306 
339.4 eon 2246 


0.47% lA | 2.96 


Clamped-Clamped Beam 





Notes: 1. All frequencies are in Hz. 


2. J # 0 indicates J assumed to eerlon mg ae 
modeling the floating element as a rectangular slab 
4 mm x 4 mm x 2 mm. 


3. The floating element mass M and (where appropriate) 
mass moment of inertia J are modeled as being lumped 
at the unclamped end of the beam element. 


4. Clamped-clamped data and verification of clamped-pinned 
calculations are from Blevins (1979). 





APPENDIX C 


UNCERTAINTY ANALYSIS 


From equation 2.1, 


8 Ke dq 
T= ——— (V- V_). 
WM ane GC O (C.1) 
e oO 
It is shown in Cook and Rabinowicz (1963) that if 
q = Fla,) (Ga2) 


where F is a function of n independent variables a, , then the 


most probable mean square error “ in the quantity q is 


S) IF \2 
~~? 2 


Sie 





assuming the a, are statistically independent. Each term 


er denotes the error or uncertainty in the quantity ¢. From 


equation C.3, 





2 
Me 
2 2 2 2 2. ames 
~ + €. + €o + €2 FT € € 

= ae d f f o\t Ss. (C4) 
e. 

where ae and the uncertainty in V(t) has been neglected since 
i 


it is an order of magnitude smaller than om 


3wt 2 
Since k = 3EI/L° ecm and A = we , equation C.3 obtains 





EZ 


i 2 2 2 
EF + ae > (3e,) + (3e,) 


2 


(OS eter 5) 9 (-01)° + 9(.05)2 =.0284 


and 


so that, from equation C.4, 


3a = 0.0284 + (.05)7 + 0.005 + (.01)7 + (.05)7 + 


lye ee EN 2 = 9 429 
Gipayes4 v/pad 


and 


ean 000 
i 


(C. 5a) 


(Gob) 


(C26) 


(G7) 


(C.8) 
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